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ABSTRACT 

This paper proposes a new phenomenology for strong incompressible MHD tur- 
bulence with nonzero cross helicity. This phenomenology is then developed into a 
quantitative Fokker-Planck model that describes the time evolution of the anisotropic 
power spectra of the fluctuations propagating parallel and anti-parallel to the back- 
ground magnetic field Bq. It is found that in steady state the power spectra of the 
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magnetic field and total energy are steeper than ' and become increasingly steep 

d&C /t, increases, where C = J d\ v ■ 5 is the cross helicity, £ is the fluctuation en- 
ergy, and is the wavevector component perpendicular to Bq. Increasing C with 
fixed £ increases the time required for energy to cascade to smaller scales, reduces 
the cascade power, and increases the anisotropy of the small-scale fluctuations. The 
implications of these results for the solar wind and solar corona are discussed in some 
detail. 

Subject headings: turbulence — magnetic fields — magnetohydrodynamics — solar 
wind — solar corona — solar flares 



Introduction 



Much of our current understanding of incompressible magnetohydrodynamic (MHD) turbu- 
lence has its roots in the pioneering work of Iroshnikov (1963) and Kraichnan (1965). These 
studies emphasized the important fact that Alfven waves travelling in the same direction along a 
background magnetic field do not interact with one another and explained how one can think of the 
cascade of energy to small scales as resulting from collisions between oppositely directed Alfven 
wave packets. They also argued that in the absence of a mean magnetic field, the magnetic field of 
the energy-containing eddies at scale L affects fluctuations on scales ^ L much in the same way 
as would a truly uniform mean magnetic field. 
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Another foundation of our current understanding is the finding that MHD turbulence is inher- 
ently anisotropic. Montgomery & Turner (1981) and Shebalin, Matthaeus, & Montgomery (1983) 
showed that a strong uniform mean magnetic field Bq inhibits the cascade of energy to small scales 
measured in the direction parallel to Bq. This early work was substantially elaborated upon by 
Higdon (1984), Goldreich & Sridhar (1995, 1997), Montgomery & Matthaeus (1995), Ng & Bhat- 
tacharjee (1996, 1997), Galtier et al (2000), Cho & Lazarian (2003), Oughton et al (2006), and 
many others. For example, Cho & Vishniac (2000) used numerical simulations to show that when 
the fluctuating magnetic field 55 is > Bq the small-scale turbulent eddies become elongated along 
the local magnetic field direction. Goldreich & Sridhar (1995) introduced the important and influ- 
ential idea of "critical balance," which holds that at each scale the linear wave period for the bulk 
of the fluctuation energy is comparable to the time for the fluctuation energy to cascade to smaller 
scales. Goldreich & Sridhar (1995, 1997), Maron & Goldreich (2001), and Lithwick & Goldre- 
ich (2001) clarified a number of important physical processes in anisotropic MHD turbulence and 
used the concept of critical balance to determine the ratio of the dimensions of turbulent eddies in 
the directions parallel and perpendicular to the local magnetic field. 

Over the last several years, research on MHD turbulence has been proceeding along several 
different lines. For example, one group of studies has attempted to determine the power spectrum, 
intermittency, and anisotropy of strong incompressible MHD turbulence using direct numerical 
simulations. (See, e.g., Cho & Vishniac 2000, MuUer & Biskamp 2000, Maron & Goldreich 2001, 
Cho et al 2002, Haugen et al 2004, MuUer & Grappin 2005, Mininni & Pouquet 2007, Perez & 
Boldyrev 2008). Another series of papers has explored the properties of anisotropic turbulence 
in weakly coUisional magnetized plasmas using gyrokinetics, a low-frequency expansion of the 
Vlasov equation that averages over the gyromotion of the particles. (Howes et al 2006, 2007a, 
2007b; Schekochihin et al 2007). These authors investigated the transition between the Alfven- 
wave cascade and a kinetic-Alfven-wave cascade at length scales of order the proton gyroradius p,, 
as well as the physics of energy dissipation and entropy production in the low-coUisionality regime. 
Turbulence at scales < p, has also been examined both numerically and analytically within the 
framework of fluid models, in particular Hall MHD and electron MHD. (Biskamp, Schwarz, & 
Drake 1996, Biskamp et al 1999, Matthaeus et al 2003; Galtier & Bhattacharjee 2003, 2005; 
Cho & Lazarian 2004; Brodin et al 2006, Shukla et al 2006). Another group of studies has in- 
vestigated the power spectrum, intermittency, and decay time of compressible MHD turbulence. 
(Oughton et al 1995, Stone et al 1998, Lithwick & Goldreich 2001, Boldyrev et al 2002, Padoan 
et al 2004, Elmegreen & Scalo 2004). Additional work by Kuznetsov (2001), Cho & Lazarian 
(2002, 2003), Chandran (2005), and Luo & Melrose (2006) has begun to address the way in which 
Alfven waves, fast magnetosonic waves, and slow magnetosonic waves interact in compressible 
weak MHD turbulence. Another recent development is the finding that strong incompressible 
MHD turbulence leads to alternating patches of alignment and anti-alignment between the veloc- 
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ity and magnetic-field fluctuations. (Boldyrev 2005, 2006; Beresnyak & Lazarian 2006; Mason, 
Cattaneo, & Boldyrev 2006; Matthaeus et al 2007) These studies examined how the degree of local 
alignment (and anti- alignment) depends upon length scale, as well as the effects of alignment upon 
the energy cascade time and the power spectrum of the turbulence. 

The topic addressed in this paper is the role of cross helicity in incompressible MHD turbu- 
lence. The cross helicity is defined as 



where v is the velocity and B is the magnetic field. The cross helicity is conserved in the absence 
of dissipation and can be thought of as a measure of the linkages between lines of vorticity and 
magnetic field lines, both of which are frozen to the fluid flow in the absence of dissipation (Moffatt 
1978). In the presence of a background magnetic field, Bq = Bqz, the cross helicity is also a 
measure of the difference between the energy of fluctuations travelling in the —z and +z directions. 
Dobrowolny, Mangeney, & Veltri (1980) showed that MHD turbulence with cross helicity decays 
to a maximally aligned state, with 5v = ±55/-\/47rp, where 5v and 55 are the fluctuating velocity 
and magnetic field and p is the mass density. Different decay rates for the energy and cross helicity 
were also demonstrated by Matthaeus & Montgomery (1980). In another early study, Grappin, 
Pouquet, & Leorat (1983) used a statistical closure, the eddy-damped quasi-normal Markovian 
(EDQNM) approximation, to study strong 3D incompressible MHD turbulence with cross helicity, 
assuming isotropic power spectra. They found that when C 7^ 0, the total energy spectrum is steeper 
than the isotropic Iroshnikov-Kraichnan k^^^^ spectrum. Pouquet et al (1988) found a similar result 
in direct numerical simulations of 2D incompressible MHD turbulence. More recently, Lithwick, 
Goldreich, & Sridhar (2007) and Beresnyak & Lazarian (2007) addressed the role of cross helicity 
in strong MHD turbulence taking into account the effects of anisotropy. 

This paper presents a new phenomenology for strong, anisotropic, incompressible MHD tur- 
bulence with nonzero cross helicity, and is organized as follows. Section [2]presents some relevant 
theoretical background. Section |3] introduces the new phenomenology as well as two nonlinear 
advection-diffusion equations that model the time evolution of the power spectra. Analytic and 
numerical solutions to this equation in the weak-turbulence and strong-turbulence regimes are pre- 
sented in Sections |4] and |5l Section \5\ also presents a simple phenomenological derivation of the 
power spectra and anisotropy in strong MHD turbulence. Section [6] presents a numerical solution 
to the advection-diffusion equation that shows the smooth transition between the weak and strong 
turbulence regimes. Section|7]addresses the case in which the parallel correlation lengths of waves 
propagating in opposite directions along the background magnetic field are unequal at the outer 
scale. In Section [8l the proposed phenomenology is applied to turbulence in the solar wind and so- 
lar corona, and in Section |9]the results of this work are compared to the recent studies of Lithwick, 
Goldreich, & Sridhar (2007) and Beresnyak & Lazarian (2007). 




(1) 
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2. Energy Cascade from Wave-Packet Collisions 



The equations of ideal incompressible MHD can be written 



+ (w^TVAi)-vw^ = -vn 



(2) 



where w = v ± (dB/y/A-np) are the Elsasser variables, v is the fluid velocity, 55 is the magnetic 
field fluctuation, p is the mass density, which is taken to be uniform and constant, va= Bq/ ^^4np 
is the Alfven speed, Bq = Bqz is the mean magnetic field, and n = (/> + 5^/87r)/p, which is 
determined by the incompressibility condition, V ■ = 0. Throughout this paper it is assumed 
that 55 < Bq and <va- 

In the limit of small- amplitude fluctuations (w"^ -C va), the nonlinear term ■ Vw^ in equa- 
tion ^ can be neglected to a first approximation, and the curl of equation Q becomes 



which is solved by setting V x w equal to an arbitrary function of z±VAt. Thus, w represents 
fluctuations with v = ±b that propagate in =i=z direction at speed va in the absence of nonlinear 
interactions. In the absence of an average velocity, the cross helicity defined in equation ([T]) can be 
rewritten as 



The cross helicity is thus proportional to the difference in energy between fluctuations propagating 
in the —z and +z directions. 

Equation Q shows that the nonlinear term is nonzero only at those locations where both w+ 
and are nonzero. Nonlinear interactions can thus be thought of as collisions between oppo- 
sitely directed wave packets (Kraichnan 1965). When both and are nonzero, equation ^ 
indicates that the fluctuations are advected not at the uniform velocity Tvaz, but rather at 
the non-uniform velocity ^vaz + w^. Maron & Goldreich (2001) elaborated upon this idea by 
showing that to lowest order in fluctuation amplitude, if one neglects the pressure term, then w+ 
wave packets are advected along the hypothetical magnetic field lines corresponding to the sum 
of Bq and the part of 55 arising from the w fluctuations. This result can be used to construct a 
geometrical picture for how wave-packet collisions cause energy to cascade to smaller scales, as 
depicted in Figure[IJ In this figure, two oppositely directed wave packets of dimension ~ X± in the 
plane perpendicular to Bq and length along Bq pass through one another and get sheared. Col- 
lisions between wavepackets of similar X± are usually the dominant mechanism for transferring 
energy from large scales to small scales. The duration of the collision illustrated in the figure is 
approximately At ~ X[|/va. The fluctuating velocity and magnetic field are taken to be in the plane 
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(4) 
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perpendicular to Bq, as is the case for linear shear Alfven waves. The magnitude of the nonlinear 

,± 



term in equation ^ is then ~ fkj 



where is the rms amplitude of the wave packet. 



The fractional change in the v and b fields of the wave packet induced by the collision is then 
roughly 

'At 



X 



(5) 
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Fig. 1.- 

sheared. 



When two wave packets collide, each wave packet follows the field lines of the other wave packet and gets 



If this fractional change is ^ 1 for both w+ and fluctuations then neither wave packet 
is altered significantly by a single collision, and the turbulence is weak. Wave packets travel a 
distance ^ before being significantly distorted, and the fluctuations can thus be viewed as 
linear waves that are only weakly perturbed by nonlinear interactions with other waves. In the 
wave-packet collision depicted in Figure [H the right-hand side of the wave packet is altered by 
the collision in almost the same way as the left-hand side, since both sides encounter essentially the 



same w+ wave packet, since the w+ packet is changed only slightly during the collision. Changes 
to the profile of a wave packet along the magnetic field are thus weaker than changes in the profile 
of a wave packet in the plane perpendicular to Bq (Shebalin et al 1983, Ng & Bhattacharjee 1997, 
Goldreich & Sridhar 1997, Bhattacharjee & Ng 2001, Perez & Boldyrev 2008). As a resuh, in the 
weak-turbulence limit, the cascade of energy to small is much less efficient than the cascade of 
energy to small (Galtier et al 2000). 

On the other hand, if the fractional change in equation ([5]) is of order unity then a wave 
packet is distorted substantially during a single collision, and the turbulence is said to be "strong." 
In the case that the fractional change in equation Q is ~ 1 for one fluctuation type, (e.g., w^) but 
<^ 1 for the other (w+), the turbulence is still referred to as strong. It should be noted that strong 
turbulence can arise when <^ va, provided that ^ X\\ . In strong turbulence energy cascades 
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to smaller to a greater extent than in weak turbulence, but the primary direction of energy flow 
in fc-space is still to larger k±, as discussed in the next section. 



3. Anisotropic MHD Turbulence with Cross Helicity 

In order to develop an analytical model, it is convenient to work in terms of the Fourier 
transforms of the fluctuating fields, given by 



(6) 



(271) 

The three-dimensional power spectrum (k) is defined by the equation 

{w^{k)-w^{ki))=A^{k)d{k + ki), (7) 

where (. . .) denotes an ensemble average. Cylindrical symmetry about Bq is assumed, so that 
A^{k) = A^{k±,k\\), where kj_ and k\\ are the components of k perpendicular and parallel to Bq. 
The mean-square velocity associated with fluctuations is then 



(8) 



It is assumed that at each value of k± there is a parallel wave number {k±) such that (1) 
the bulk of the fluctuation energy is at \k\\ \ < fcy {k±^) and (2) A^{kj_,k\\) depends only weakly 
on for I < k^{k^). A w"^ wavepacket at perpendicular scale kj^ then has a correlation length 

in the direction of the mean field of ~ (^k^^ . The rms amplitude of the fluctuating Elsasser 
fields at a perpendicular scale denoted wf^, is given by 

{wf^)^r^A^{k^,0)klk^. (9) 

As described in the section [2l when a wave packet at scale kj^ collides with a w"^ wave packet 
at scale the fractional change in the packet resulting from the collision is approximately 

The wave number kf is defined to be the value of ^ |^ for which = 1- Thus, 
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3.1. The energy cascade time 

When ^ fc^ , the value of %^ is ^ 1 and a w+ is only weakly affected by a single collision 

with a wave packet. Each such collision requires a time (^|| v^)^^. The effects of successive 
collisions add incoherently, and thus (x^^)^^ collisions are required for the w+ wave packet to be 
strongly distorted, and for its energy to pass to smaller scales. The cascade time for a w+ wave 
packet at perpendicular scale kj^ is thus roughly 

'^t ~ (^wVAy\')CkJ^^ ~ 7— (weak turbulence). (12) 

kc Va 

Similarly, if < 1, then x^^ ~ {^VAy^. 

When ~ k^, the value of is ~ 1, a is strongly distorted during a single wave packet 
collision, and the turbulence is strong. Each such collision takes a time va)^^- Since fcy ~ k^, 

x^ ~ (^11 va)^^ ~ — (strong turbulence). (13) 

^ kc VA 

Similarly, if ~ 1, then x^^ ~ (^+va)"^ 

The case k^ <^ kf (i.e., %^ ^ 1) is explicitly excluded from the discussion. Initial conditions 

could in principle be set up in which k\^ <^ kf. However, the cascade mechanisms described in 

section 13.21 will not produce the condition ^ kf if it is not initially present. It should be 
emphasized that in both weak turbulence and strong turbulence, the cascade time is given by the 
same formula, x^ ~ (fc^VA)^^ which involves the spectrum evaluated only at k\\ = 0. 



3.2. The Cascade of Energy to Larger | 

The two basic mechanisms for transferring fluctuation energy to larger ^|| were identified 
by Lithwick, Goldreich, & Sridhar (2007). The first of these can be called "propagation with 
distortion." Suppose a w+ wave packet of perpendicular scale kj^ and arbitrarily large initial 
parallel correlation length begins colliding at ? = with a stream of wave packets of similar 
perpendicular scale. At time t = x^ , the leading edge of the w+ wave packet has been distorted 
substantially by the stream of wave packets, but the trailing portion of the w+ wave packet 
at distances > 2vaX^ behind the leading edge has not yet encountered the stream of wave 
packets. If the parallel correlation length of the w+ wave packet is initially > 2vaX^, then during 
a time x^ the w+ wave packet acquires a spatial variation in the direction of the background 
magnetic field of length scale ~ 2vaX^ ~ 2(k^)^^. This process is modeled as diffusion of 
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fluctuation energy in the k\\ direction with diffusion coefficient ~ (Afc||)^/A?, where A^y = 
and A? = x^. "Propagation with distortion" then leads to a value of of ~ {kJ^'vA- 

The second mechanism identified by Lithwick, Goldreich, & Sridhar (2007) can be called 
"uncorrelated cascade." Consider a w+ wave packet of perpendicular scale k'^ and arbitrarily 
large parallel correlation length, and consider two points within the wave packet, P\ and P2, that 
move with the wave packet at velocity —vaz and are separated by a distance along Bq of 2va\^ ~ 
2{k^)^^ . The wave packets at perpendicular scale k'^ encountered by the portions of the 
w+ wave packet at Pi and P2 are then uncorrelated, because wave packets are substantially 
distorted while propagating between Pi and P2. Thus, the way in which the wave packet 
cascades at location P\ is not correlated with the way in which the w+ wave packet cascades at 
location P2. If the parallel correlation length of the w+ wave packet is initially > 2vaX^, then 
wave-packet collisions introduce a spatial variation along Bq into the w+ wave packet of length 
scale ~ 2va\^ ~ 2(^+)^^ during a time x^. Again, we model this as diffusion of fluctuation 
energy in the fcy direction with D"^ ~ (Afc[|)^/Af and /St = x^, but now Afcy = kf. "Uncorrelated 
cascade" thus leads to a fcy -diffusion coefficient of ~ {kf)^kfvA. 

Accounting for both mechanisms, one can write 

Df ~ {kc,m.x)^k^VA, (14) 

where ^f,max(^±) is the larger of k^{k±) and k^{k±). If k^ > k^, then w+ energy diffuses in ^[| 
primarily through the "uncorrelated cascade" mechanism, while energy diffuses in ^|| primarily 
through the "propagation with distortion cascade" mechanism. 



3.3. Advection-Diffusion Model for the Power Spectra 

The phenomenology described in the preceding sections is encapsulated by the following 
nonlinear advection-diffusion equation, 

where A^ is shorthand for A^(^^,fc||), c\ and C2 are dimensionless constants of order unity, and 
5*^ and — T^A^ are forcing and damping terms, respectively. The first term on the right-hand side 
of equation (fT5l) represents advection of fluctuation energy to larger k^, while the second term 
represents diffusion of fluctuation energy to larger \k^. The quantity x^^ is an effective cascade 
time at perpendicular scale Usually, the transfer of energy to small scales is dominated by 
local interactions in ^-space, and the cascade time for a w+ wave packet is ~ (fc^v^)^'. In some 
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cases, however, the shearing of small-scale wave packets by much larger-scale wave packets can 
become important. To account for such cases, the effective cascade time is taken to be 



v-1 



max 



for < q^<k^ 



(16) 



i.e., (xir , ) ^ is the maximum value of k^VA for all perpendicular wave numbers between zero 



and k^. The flux of energy to larger is 



The term is given by 



(17) 



(18) 



and is included so that e"*^ increases as the spectrum becomes a more steeply declining function 
of k^, in accordance with weak turbulence theory (Galtier et al 2000, Lithwick & Goldreich 2003). 
To match the energy flux in weak turbulence theory in the limit of zero cross helicity, one must 
se£l 

c,=-f. (19) 



where 
7 = 



1 J-i [x^-y^p 



For simplicity. 



C2=l. 



(21) 



4. Steady-State Weak Turbulence 

This section addresses weak turbulence in which ^ii*^ ~ at the outer scale. The weak- 
turbulence condition, ^ 1, is equivalent to the condition kf <^ k^. Because /:[| -diffusion 



'The value of ci in equation (fT9] l is a factor of 2 larger than the value that follows from the results of Galtier 
et al (2000). It appears that this discrepancy results from the omission of a factor of 2 in equation (54) of Galtier 
et al (2000). This can be seen by starting from equation (46) of Galtier et al (2000) and using the expression on 
page 1045 of Leith & Kraichnan (1972) to simplify polar integrals of the form / d^pd^qd{k~ p — q)F{k,p,q) for 
two-dimensional wave vectors k, p, and q, where F is a function only of the wave-vector magnitudes and the integral 
is over all values of p and q. 
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involves a Afcy ~ kf during a time x^, the fcy -increment over which energy diffuses while cascading 
to larger is much less than the breadth of the spectrum in the k\\ direction (~ kt), so the ^ii- 



diffusion terms can be ignored to a good approximation. In this case, equation (1151) possesses a 
steady-state solution in which e+ and £^ are constant, and in which 

A^=g^{k\\)kl"\ (22) 

where g"*" and g~ are arbitrary functions of k\\ , and where 

n+-^n"==6, (23) 

with 2 < < 4. Equations (f22)) and (l23l) match the results of weak turbulence theory for in- 
compressible MHD turbulence if one allows only for three-wave interactions among shear Alfven 
waves (Galtier & Chandran 2006), or if one considers only the limit that k± ^ k\\ (Galtier et al 
2002). If one writes = 3±a with |a| < 1 and sets g^{k]^\) = g^(fc||), then equation (flTl) gives 

^ = |i^. (24) 

E- 2-a 

In the limit a ^ 1, £+/£ = 1 + a, in agreement with the weak-turbulence-theory result for 
k± ^ (Lithwick & Goldreich 2003), as in the weak-turbulence advection-diffusion model 
of Lithwick & Goldreich (2003). 

In steady state, A+(fcj^,0) and A^(A:^,0) are forced to be equal at the dissipation scale so that 
= . This phenomenon of "pinning" was discovered by Grappin et al (1983) for strong MHD 
turbulence, and further elaborated upon by Lithwick & Goldreich (2003) for the case of weak 
turbulence. The dominant fluctuation type then has the steeper spectrum. If £+/£ is fixed, then 
the ratio / of the rms amplitudes of the two fluctuation types at the outer scale k^ ^ increases 
as k^/kf increases, where k^j is the dissipation wave number. Alternatively, if >v^/w^ is fixed, 
then e+/e^ approaches unity as k^/kf oo. 

Several of these results are illustrated by the numerical solution to equation (fTSl ) shown in 
Figure[2l This solution is obtained using a logarithmic grid for k^, with k±^i — kol'^" for <i <N. 
Similarly, k^^j = kQ^/" for 1 < j < M, but k^^j = for j = 0. Af is advanced forward in time 
using a semi-implicit algorithm, in which the terms hf, x^^^ , and kf on the right-hand side of 
equation (fT5l) are evaluated at the beginning of the time step, and the terms on the right-hand 
side of equation (fT5l) are evaluated at the end of the time step. The algorithm employs operator 
splitting, treating the ^^-advection, forcing, and damping in one stage, and the -diffusion is a 
second stage. In this approach, the matrix that has to be inverted to execute each semi-implicit 
time step is tri-diagonal. An advantage of this procedure over a fully explicit method is that the 
time step is not limited by the fcii -diffusion time at large k± and small kn . The discretized equations 
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are written in terms of the energy fluxes between neighboring cells, so that in the absence of forcing 
and dissipation the algorithm conserves fluctuation energy to machine accuracy. For the numerical 
solution plotted in Figured = 80, M = 16, n = 4, = S^k^expi-k^ /kj), = 1.25^, kf = 
5kQ, and 7^ = 2k^v, where v is an effective viscosity. The initial spectra are set equal to zero, 
and the equations are integrated forward in time until a steady state is reached. In steady state, 
5v+ = 2.5 X lO^^^A and Sv" = 6.4 x IO^Sa. 




Fig. 2. — Numerical solution of equation ( fTSI ) in the weak-turbulence limit. Left panel: The dimensionless ID power 
spectrum defined in equation dZSl l. Middle panel: The spectral slopes at k\\ = 0. Right panel: The weighted value of ^|| 
defined in equation ( l26b . In all panels, the solid lines refer to w+ and the dashed lines refer to w^. In the right-hand 
panel, the two lines are almost on top of each other 

The left-hand panel of Figure [21 is a plot of the dimensionless one-dimensional power spec- 
trum, 

E^{k^) = ^ r dk\\A^{k^,k\\), (25) 

which is proportional to the energy per unit k^ in fluctuations. The middle panel of Figure [2] 
shows that in the inertial range, dlnA^ {k±,0) /dink ^ —3.2 and dlnA^ {k±,0) / dink ~ —2.8, as 
expected for — 1.2. The right-hand panel shows that the weighted value of k^^, 

/oo 
dk^i \ki\\A^{k^,k\\) 

{K^} = , (26) 

/ dk\\A^{k^,k\\) 



is roughly constant in the inertial range. 
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5. Steady-State Strong Turbulence 

This section addresses strong turbulence with %^ ~ 1, %^ < 1, > w^, and ~ fcy at 

the outer scale kj^. The discussion allows for the possibility that %^ ^ 1. As fluctuation energy 

cascades to larger k±, it diffuses to larger \k\\\, so that k^ and fcy increase with increasing k±. 
Moreover, for both w+ and w^, the fluctuation energy diffuses over a -increment of ~ k^ during 
one cascade time. For the steady-state solutions of interest, k^ is an increasing function of k±, and 
thus at each k^ we will have that k^ ~ fcy ~ k^. One can thus define a single parallel-wavenumber, 
{k±), to describe the spectra, with 

~ ^ii'' ~ ~ k^ (27) 

at each k^. Since k^ ~ k^ at each scale, 

~ 1 (28) 

throughout the inertial range. On the other hand, since can be much less than w^, %^ can 
be < 1. 

The cascade time for the w~ fluctuations is given by the strong-turbulence phenomenology of 
equation (fT3l) . so that the energy flux in fluctuations is 

e- r^thL^k^w+ (w,)2_ (29) 

Allowing for the possibility that %^ <^ 1, the cascade time of the fluctuations follows the 
weak- turbulence phenomenology of equation (fT2l) . This formula is also accurate for as large 
as 1 (in which case ~ ^^^)- The energy flux in w+ fluctuations is then 

£+ ~ y ~ (W. )2, (30) 

which is roughly the same as e . It is assumed that the energy flux depends on the spectral slope 
as in weak turbulence, so that the fluctuation type with the steeper spectrum has the larger energy 
flux. If 

w^ock-^"^^ (31) 
then equations (|29l) and (|30l) imply that when £+ and £^ are independent of k^, 

a+ + 2a" = l. (32) 
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The condition that ~ 1 throughout the inertial range then implies that 

^ock^-"^. (33) 

As discussed in earlier studies (Grappin et al 1983, Lithwick & Goldreich 2003), the spectra are 
pinned at the dissipation scale, so that the dominant fluctuation type will have the steeper spectrum 
and a somewhat larger energy flux. For the zero-cross-helicity case, equations (|32l ) and (|33] ) give 
~ °^ k^^^^ and oc k^^, as in the work of Goldreich & Sridhar (1995) [see also Higdon 
(1984)]. ^ 

When =yf = 0, equation ([T?] ) possesses an analytical solution that reproduces the above 
scalings. This solution can be obtained by starting with the assumptions that 

A^{k^,0) = cfkl''\ (34) 

that the energy cascade is dominated by local interactions, and that A+(fcj^,0) > A^{kj_,0) for 
all k^. Equation (fT6l ) then becomes (x^^^ )^^ = k^A^{k^,0) /va, and fcc,max = k^- Upon defining 



and 



one can rewrite equation (fT5]) as 



with 



= k'-'^Af (35) 



5 = yt^-2fc+^ (3g) 



Jk _j^±Jk^ (37) 



ds dkl 



= .riw^ ..4 - (38) 



C2(c+)2 

ci{%-2b+){b^-l)v\ 



Equation (|37l) is solved by taking 



Requiring that equation (|34l ) be satisfied, one finds that = cf and 

2b+ + b- = 10. (40) 

The dominance of local interactions requires that b^ < 4, and thus > 2. When forcing and 
dissipation are taken into account, the exact solution becomes an approximate solution that is valid 
only within the inertial range. In this case, > because the spectra are pinned at the dissipation 
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scale whereas A+(fc^,0) is larger than A (fc^,0) within the inertial range. Equation (l39l) implies 
that 

_, __ c+k^-^^ 

r^VWs^^^ — , (41) 

II II ^2 

where the dimensionless constants in the expression for have been dropped, but c^, which has 
dimensions, has been kept. Equations (fTTI) . (I34l) . and (|4T]) show that fc+ ~ for all fc^, so that 
~ 1 for all k^. Equation (|25] ) gives 

oc A±(A:i,0)/ti^if oc A:^-^^-**, (42) 

from which it follows that 



and 



<-^i-'", (43) 

A6-b+-b-)/2 

w^^ock\ " . (44) 

This solution reduces to the critical-balance solution of Goldreich & Sridhar (1995) when = 

-[- 2/3 -I- 1/3 

= 10/3, in which case fcy k^ , and oc ' . Comparing equations (l43l) and (|44l) with 
equation (|3T| ). it can be seen that b^ corresponds to 3 + a+ and b corresponds to 3 + 2a^ — a~^. 
Equation (|4TI) is thus equivalent to equation (l33l) . and equation (|40|) is equivalent to equation (|32|) . 

Figure |3] shows the results from a numerical solution of equation ([T?] ), obtained by integrating 
forward in time as described in section |4] with the spectra initially equal to zero. The numerical so- 
lution was obtained by setting = Sq k^ exp{—k^ / kj:) with Sq = 1.25q and using the parameters 
(defined in section |4l) N = SO, M = 60, n = 4, kf = 5ko, and = Ik^v, where the constant v is an 
effective viscosity. The rms velocities at steady state are 5v+ = 1.5va and 5v^ = O.IOva. The left- 
hand panel shows that the one-dimensional energy spectrum E'^ is k^^'^^ in the inertial range, 
which corresponds to b^ = 3.57 in equation (|42l) . Equation (|40l ) then gives b = 2.86. The dotted 
lines in the middle panel of Figure |3] correspond to the values of b^ = 3.57 and b^ = 2.86, which 
are reasonably close to the values of — JlnA'*^(^^,0)/Jln^^ in the numerical solution, although 
these latter values vary throughout the inertial range in the numerical solution. For b^ = 3.57, 
equation (|4T]) gives k\^ k^^^, which is a close match to the numerical solution, as shown in the 
right-hand panel of Figure [3l The left-hand panel of Figure [3] shows that the steady-state solutions 
for A+ and are "pinned" at the dissipation scale, as expected. 

It should be noted that when ~ 1 and <^ 1, the dominant w+ fluctuations are only 
weakly damped by nonlinear interactions with waves, in the sense that is much larger 
than the linear wave period. On the other hand, for the smaller-amplitude w fluctuations, the 
linear wave period and cascade time are comparable. Thus, paradoxically, the larger-amplitude w+ 
fluctuations can be described as waves, or, more precisely, a non-sinusoidal wave train, whereas 
the smaller-amplitude fluctuations can not be accurately described as waves. 
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Fig. 3. — Numerical solution of equation ( fTsT i for strong turbulence with ~ 1- Left panel: The dimensionless 
ID power spectrum defined in equation (IZST l. Middle panel: The spectral slopes at k\\ ~ 0. Right panel: The weighted 
value of defined in equation (|26] |. In all panels, the solid lines refer to w+ and the dashed lines refer to . 

6. Transition Between Weak Turbulence and Strong Turbulence 

This section again addresses turbulence in which ~ fcy at the outer-scale wavenumber, kf. 
In the weak-turbulence limit, and increase with increasing k^. If the dissipation wavenum- 
ber kd is sufficiently large, then %^ and/or will increase to a value of order unity at some k^ 
within the inertial range. This perpendicular wavenumber is denoted ^trans- The turbulence will 
then be described by the weak- turbulence scalings of section |4] for kf <^ k± <^ ^trans^ and by the 
strong-turbulence scalings of section [5] for fctrans <^k± -C fc^. Figure |4] shows a numerical solution 
of equation (fTSi) that illustrates how the turbulence makes this transition in a smooth manner. At 
small wavenumbers, this solution is similar to the weak-turbulence solution plotted in Figure [21 
and at large wavenumbers it is similar to the strong-turbulence solution plotted in Figure [H The 
solution shown in Figure H] was obtained by integrating equation (fTSi) forward in time to steady 
state using the numerical method described in section |4l The spectra were initially set equal to 
zero. The numerical solution was obtained by setting 5^ = 5q fc^exp(— fc^/fc^^) with Sq = 1.2Sq 
and using the parameters = 100, M = 56, n = 4, kf = 5ko, and 7^ = 2k^v, where the constant V 
is an effective viscosity. The rms velocities at steady state are 5v+ = 0.32va and 6v^ = 0.012va. 



7. Unequal Parallel Correlation Lengths at the Outer Scale 

In sections m through [6l it was assumed that k\^ ~ fcy at the outer scale. This assumption is 
applicable to many settings. For example, in a plasma of dimension L that is stirred by a force 
that has a correlation length I <^L, the velocity fluctuations that are excited have a correlation 
length /, and this correlation length is imprinted on both the w+ and fluctuations. On the other 
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Fig. 4. — Numerical solution of equation ( fTSl l showing a smooth transition from the weak-turbulence limit at 
small k± to the strong-turbulence limit at large k±. Left panel: The dimensionless ID power spectrum defined in 
equation dZSl l. Middle panel: The spectral slopes at = 0. Right panel: The weighted value of ^|| defined in equa- 
tion (|26] |. In all panels, the solid lines refer to and the dashed lines refer to w^. 

hand, if waves are launched along the magnetic field into a bounded plasma from opposite sides 
of the plasma, and the waves from one side have a much larger parallel correlation length than the 
waves from the other side, it is possible to set up turbulence in which the two wave types have very 
different parallel correlation lengths at the outer scale. This situation is discussed briefly in this 
section. 

For strong turbulence, if both %^ and %^ are ~ 1 at some perpendicular scale , but one 
fluctuation type, say w+, has a much smaller parallel correlation length than the other (and thus 
a much larger amplitude), then during a time x^^ the "propagation with distortion" mechanism 

discussed in section [3^ will increase fcy until it equals , which will cause %^ to become ^ 1 at 
scale k'^ . At smaller scales, the solution can be described by the scalings presented in section [51 
in which k^{k^) ~ {k^). Similarly, if %^ ~ 1 but ^ 1 at some scale fcj^ and if ^||'^ ^ 
at that scale, then during a time the "propagation with distortion" mechanism discussed in 

section [X2] will again increase k\^ until it equals ^||^, the parallel scales will remain comparable at 
smaller perpendicular scales, and the solution can be described by the scalings in section [51 The 
case in which ~ 1, -C 1, and k'^ ^ fcy is not addressed in this paper. 



8. Implications for Turbulence in the Solar Corona and Solar Wind 

In this section, the preceding analysis of incompressible MHD turbulence is applied to the 
solar wind and solar corona. It should be noted at the outset, however, that the solar wind and 
solar corona (beyond roughly r = 1.5i?0, where r is distance from the Sun's center) are in the 
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coUisionless regime, and the pressure tensor is not isotropic as assumed in ideal MHD. Moreover, 
the value of p = 87rj!?/5^ is <^ 1 in the corona and typically ~ 1 in the solar wind at 1 AU, whereas 
incompressible MHD corresponds to the limit |3 — ^ oo. A preliminary question that needs to be 
addressed is thus the extent to which incompressible MHD is an accurate model for these plasmas. 

Schekochihin et al (2007) have recently carried out extensive calculations based on kinetic 
theory that provide a detailed answer to this question. These authors examined anisotropic turbu- 
lence in weakly coUisional magnetized plasmas using gyrokinetics, a low-frequency expansion of 
the Vlasov equation that averages over the gyromotion of the particles. By applying the form of 
the gyrokinetic expansion derived by Howes et al (2006), Schekochihin et al (2007) showed ana- 
lytically that non-compressive Alfvenic turbulence in the quasi-2D regime (i.e., k± ^) can be 
accurately described using reduced MHD in both the coUisional and coUisionless limits, regardless 
of p, provided that the length scales of the fluctuations are much larger than the proton gyroradius 
and the frequencies are much less than the proton cyclotron frequency. Since non-compressive 
quasi-2D fluctuations are thought to be the dominant component of the turbulence in the solar 
wind (see, e.g., Bieber et al 1994) and the solar corona (Dmitruk & Matthaeus 2003, Cranmer & 
van Ballegooijen 2005), incompressible MHD is a useful approximation for modeling turbulence 
in these settings. 



8.1. Cross helicity in the solar wind and solar corona 

Cross helicity in the solar wind has been measured in situ by several different spacecraft. In 
terms of the Elsasser variables w"^, there is a substantial excess of outward propagating fluctuations 
(taken to be w+ throughout this section) over inward propagating fluctuations (taken to be w~) in 
the inner heliosphere, although this imbalance decreases with increasing r, as seen, for example, 
in Voyager data for low heliographic latitude (Matthaeus & Goldstein 1982, Roberts et al 1987) 
and Ulysses data at high latitude (Goldstein et al 1995). In a study of Ulysses and Helios data, 
Bavassano et al (2000) found that oc r'^-^'^ for r < 2.6 AU and e+/e~ ~ 2 for 3 AU < r < 

5 AU, where is the energy per unit mass associated with fluctuations. These numbers are 
intended as illustrative average values, as individual measurements of in the solar wind vary 
significantly. 

Although it has not been directly measured, the ratio e'^/e~ is likely very large in open-field- 
line regions of the solar corona. This can be seen from the work of Cranmer & van Ballegooijen 
(2005), who modeled the generation of Alfven waves by the observed motions of field-line foot- 
points in the photosphere, and the propagation and reflection of these waves as they travel along 
open field lines from the photosphere out into the interplanetary medium. They found that the 
ratio of the frequency-integrated rms Elsasser variables (w"*" and w~) is ~ 30 at r = 2Rq (i.e.. 
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~ 900). Verdini & Velli (2007) developed a different model for the propagation, reflection, 
and turbulent dissipation of Alfven waves in the solar atmosphere and solar wind and found that 
~ 80 at r = 2Rq. Based on these results, one can make the rough estimate that 

^-10 (atr = 2i?0), (45) 

where kf is the perpendicular wavenumber at the outer scale. 



8.2. Is quasi-2D turbulence in the corona and solar wind weak or strong? 

In much of the solar wind, 55 is comparable to and the turbulence is in the strong- 
turbulence regime with x^, ~ 1. For the corona, Cranmer & van Ballegooijen (2005) found 
that the outer-scale fluctuations in open-field-line regions have periods of T = 1 — 5 minutes, 
5v ~ 100 km/s, and perpendicular correlation lengths of ~ fcj^ ~ 10^ km. The Alfven speed in 
their model corona is between 2000 and 3000 km/s at r = 2Rq, and thus 5v ^ v^. However, the 
parallel correlation length L|| of the outer-scale fluctuations is ~ vaT = 1 .8 — 9 x 10^ km, which is 
^ L±. Because L||/L^ ~ va/6v, 

it, - 1 (46) 

and the low-frequency fluctuations launched into the corona by footpoint motions are in the strong- 
turbulence regime. There may be an additional population of higher- frequency waves in the weak- 
turbulence regime, but these are not discussed here. 



8.3. Parallel correlation lengths of inward and outward waves 

In open-field-line regions of the corona, when w+ waves are reflected, the resulting w~ waves 
have the same frequencies as the w+ waves. On the other hand, wave-reflection is more efficient 
at lower frequencies (Velli 1993), so if there is a range of wave frequencies at each k±, the energy- 
weighted average frequency of inward waves would tend to be somewhat lower than that of the 
outward waves. This suggests that at the outer scale the parallel correlation length L|| of the 
fluctuations is somewhat larger than the value of Ly for the w+ fluctuations. However, given 
equation (|46|) . the w+ fluctuations imprint their parallel correlation length on the w fluctuations 
during a single turnover time x^, as argued in section |7] The parallel correlation lengths of the w+ 
and w fluctuations in the corona can thus be taken to be approximately equal at the outer scale, 
and hence also at smaller scales. The same approximation is reasonable for turbulence in the solar 
wind. 



- 19- 



8.4. Energy Dissipation Rate 

If we take kf to he the perpendicular wave number at the outer scale, w^^ to be the rms ampli- 
tude of the outward-propagating fluctuations at the outer scale, and to be the rms amplitude of 
the Sunward-propagating fluctuations at the outer scale, then equations (|29l) , (|30l) . and (|46l ) imply 
that 

e+ ~ £- ~ kfw^j.{wj^jf. (47) 

This estimate is a factor of ~ w^^/^^^y smaller than the standard strong-turbulence estimate of 
e+ ~ kf{w^ '^hat appears in many studies (e.g., Zhou & Matthaeus 1990, Cranmer & Van 
Ballegooijen 2005, Lithwick, Goldreich, & Sridhar 2007, Verdini & Velli 2007). This difference 
has important implications for turbulent heating of the solar corona and solar wind. 



8.5. Cascade Time 

For the energetically dominant fluctuations, the cascade time is much longer than 
the linear wave period, at least at scales much larger than the dissipation scale. This result is 
important for determining the conditions under which turbulence can be a viable mechanism to 
explain the heating of the solar corona. Observations taken with the Ultraviolet Coronagraph 
Spectrometer (UVCS) indicate that there is strong heating of coronal plasma at r < 2Rq (Kohl et 
al 1998, Antonucci et al 2000). An appealing model to explain this heating is that low-frequency 
Alfven waves are launched by turbulent motions of field-line footpoints in the photosphere, that 
some of these waves are reflected, and that interactions between oppositely directed Alfven wave 
packets in the corona causes the wave energy to cascade to small scales and dissipate (Matthaeus 
et al 1999, 2002; Dmitruk et al 2001, 2002; Cranmer & van Ballegooijen 2005, 2007; Verdini 
& Velli 2007). In one version of this model, the waves that cross the transition region and enter 
the corona have not yet undergone a turbulent cascade, and their energy is concentrated at the 
fairly long periods (> 1 minute) characteristic of the observed footpoint motions that are believed 
to make the dominant contribution to the outward directed wave flux. In order for this scenario 
to explain the UVCS measurements, there needs to be time for the outer-scale waves to cascade 
within the corona before they travel beyond r ~ 2Rq. If, as above, we take fc/ to be the value of 
at the outer scale, to be the rms amplitude of the outward waves at the outer scale, and L\\ to be 
the parallel correlation length of the fluctuations at the outer scale, then equation (fT2l) can be used 
to express the cascade time for the outward waves at the outer scale as 

<r-{—] (48) 
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where %^ ~ kfW^L\^\/vA- Thus, for waves with a period L\^/va ~ 1 minute, equations (|45l) . (|46l) . 
and (l48l) give x^^ ~ 100 minutes. On the other hand, the Alfven speed in a coronal hole at r < 2Rq 
is ~ 2000 — 3000 km/s (Cranmer & van Ballegooijen 2005), and the time for an Alfven wave to 
travel from the coronal base out to r = 2Rq is 4-6 minutes. There is thus not enough time for the 
energy of waves with periods > 1 minute to cascade and dissipate within a few solar radii of the 
Sun. 

Dmitruk & Matthaeus (2003) and Verdini & Velli (2007) avoid this difficulty by postulating 
that a broad frequency spectrum of waves is launched upwards from the photosphere. Another pos- 
sible way around this difficulty is the development of a broad frequency spectrum of fluctuations 
from wave-packet collisions in the chromosphere, in which the energies of inward and outward 
propagating waves are comparable due to strong wave reflection at the transition region. 

8.6. Spectral Index 

Much of the discussion of the inertial-range power spectrum of solar-wind turbulence has 
focused on the question of whether the spectral index is closer to the Kolmogorov (1941) value 
of —5/3 or the Iroshnikov-Kraichnan value of —3/2 (h^oshnikov 1963, Kraichnan 1965). A value 
of —5/3 is supported by a number of theoretical studies (e.g., Montgomery & Turner 1981, Higdon 
1984, Goldreich & Sridhar 1995) and numerical simulations (Cho & Vishniac 2000, Miiller & 
Biskamp 2000, Cho et al 2002, Haugen et al 2004). A value value of —3/2 is supported by a 
second group of theoretical studies (Boldyrev 2005, 2006; Mason et al 2006; see also Beresnyak 
& Lazarian 2006) and numerical simulations (Maron & Goldreich 2001, Miiller et al 2003, Miiller 
& Grappin 2005, Mininni & Pouquet 2007). It should be noted that all of the above-mentioned 
studies address MHD turbulence with negligible cross helicity. 

Spacecraft measurements yield frequency spectra for the magnetic field and velocity fluctu- 
ations, where the frequency / is approximately kyU /In, where kr is the radial component of the 
wave-vector and U is the solar-wind speed [Taylor's (1938) "frozen-in flow hypothesis"]. Below a 
spectral-break frequency fi,, the spectra are typically fairly flat, being approximately proportional 
to /^^ (Matthaeus & Goldstein 1986). At / > Z^, the spectra steepen. The time scale correspond- 
ing to the spectral break, increases with increasing r. For example, Bruno & Carbone (2005) 
found that f^^ was 0.06 hours at 0.3 AU, 0.16 hours at 0.7 AU, and 0.4 hours at 0.9 AU in a sample 
of Helios 2 data. In two other studies based on data from several spacecraft, Matthaeus & Gold- 
stein (1986) found that /^"^ ~ 3.5 hours at 1 AU, while Klein et al (1992) found /^^^ ~ 12 hours at 
4 AU. The inertial range roughly corresponds to frequencies in the interval fb<f< fd^ where 
is the frequency corresponding to the dissipation scale. At 1 AU, fd ~ 0.3 s^^ (Smith et al 2006). 
A large number of inertial-range spectral indices have been reported in the literature. For example. 
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Matthaeus & Goldstein (1982) found a spectral index of —1.73 ±0.08 for the magnetic field in 
Voyager data at r = 1 AU, and a spectral index of —1.69 ±0.08 for the total energy. Goldstein 
et al (1995, Fig. 1) found that the spectral index for the w+ fluctuations was slightly steeper than 
—5/3 in Ulysses data at 2 AU and 4 AU. Their results also suggest a shallower spectrum, con- 
sistent with the idea that the spectra are pinned at the dissipation wavenumber kj. In a study of 
Helios 2 magnetic-field data, Bruno & Carbone (2005, Figure 23) found spectral indices of —1.72 
at 0.3 AU, -1.67 at 0.7 AU, and -1.70 at 0.9 AU. Using data from the WIND spacecraft at 1 AU, 
Podesta et al (2007) found a total-energy spectral index of — 1 .63, with the velocity spectrum flatter 
than the magnetic spectrum. Marsch & Tu (1996) found a spectral index of —1.65 ±0.01 for the 
magnetic field at 1 AU in Helios 2 data. Horbury & Balogh (1995) found a spectral index close 
to —5/3 for the magnetic field in Ulysses data at 2.5 AU. Using magnetic-field data from the ACE 
spacecraft at 1 AU, Smith et al (2006) found a spectral index of —1.63 ±0.14 in open-field-line 
regions and —1.56 ±0.16 in magnetic clouds. Smith (2003) found spectral indices between —1.7 
and —1.8 in a study of Ulysses magnetic-field data covering a range of heliographic latitudes and 
radii. 

Overall, the spectra are more consistent with a Kolmogorov scaling than an Iroshnikov- 
Kraichnan scaling. It should be emphasized, however, that the observations in several cases are 
consistent with inertial range spectra that are steeper than a Kolmogorov spectrum. Spectral in- 
dices > 5/3 have been found in previous theoretical studies of weak incompressible MHD turbu- 
lence (Bhattacharjee & Ng 1997, Goldreich & Sridhar 1997, Galtier et al 2000, Bhattacharjee & 
Ng 2001, Perez & Boldyrev 2008), as well as strong isotropic MHD turbulence with cross helicity 
(Grappin et al 1983). In this paper, it is argued that spectral indices > 5/3 are a consequence of 
cross helicity in strong anisotropic MHD turbulence. 

The simplest way to apply this paper to solar wind turbulence is to model the solar wind 
fluctuations at some location as steady-state, forced, homogeneous turbulence with the same av- 
erage value of w^^/w^^, where w^^ and are the rms amplitudes of the fluctuations at 

the outer scale kj^. Upon setting kj" , one can write (>^^/>v^)^ °^ {k±/kj)^^"^^^" , 

where it is assumed that the spectra are equal at the dissipation wave number k^. Equation (|32l) 
then gives (w^^/w^^)^ ~ (k^/k^y^^^^ , and the value of a+ can be obtained from the equation 

(w^/w^^.)^ — (kf/kd)^^^"^. The total-energy spectrum, E{k±) = E^{k±) +E^{k±), is approxi- 
mately A:^^(w^)^, although it is flatter than fc^^(w^)^ near the dissipation scale where the flatter 
spectrum of the w fluctuations is important. Thus, at scales much larger than the dissipation scale. 



where 



E{k^)ock^'^, (49) 
5 21ogio[(w+)VKP'] 

^ + T^ /, /, / • (50) 

3 3log^Q{kj/kf) 
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Upon defining the outer-scale fractional cross helicity as 

Or 



(51) 



one can rewrite equation (pOl) as 



^5 ^ 21ogio[(l+o,)/(l-o, 
3 3logiQ{kd/kf) 



(52) 



The spectral index from equation (|52l) is plotted in Figure[5l assuming that kj/kf = fd/ ft = 3780, 
where ft, = (3.5 hours)^^ is the break frequency at 1 AU discussed above (Matthaeus & Goldstein 
1986) and = 0.3 s^^ is the frequency at the dissipation scale. (Smith et al 2006) With this 
choice, <7 = 1.78 for {w^j)^ / {w^^)^ = 4 and <? = 1.74 for (w^^.)^/(w^^)^ = 2. When equation ^ 
is applied to the solar wind, Oc should be interpreted as the cross helicity at the outer scale kj^ 
averaged over at least a few outer-scale fluctuations. 
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Fig. 5. — Dependence of spectral index of the total energy spectrum, q, on the outer-scale fractional cross helicity Oc- 
The ratio of the dissipation wavenumber kj to the perpendicular wavenumber at the outer-scale kf in equation ( |52] | is 
taken to be 3780. 



Some caution, however, is warranted when applying equation (1521) to the solar wind because 
solar wind conditions vary with r, while equation (|52l) is based on results for homogeneous turbu- 
lence. The condition under which it is valid to treat the solar- wind fluctuations at a some r with an 
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observed Oc as homogeneous turbulence with the same Oc is that the cascade time at the outer scale 
be much smaller than the time ?adv — r/U for turbulent structures to be advected a distance r, 
where U is the solar wind speed. However, this condition is often either violated or only marginally 
satisfied. This is illustrated by the results of Matthaeus & Goldstein (1982) based on four days of 
Voyager data (in their Table 1). For the fluctuations at 1 AU, these authors found an rms velocity 
fluctuation of 5v = 25.5 km/s [which is comparable to va at 1 AU], a velocity correlation length 
(outer scale) of Lj^ ~ 2.83 x 10^^ cm, and an average solar-wind speed of 352 km/s. In this par- 
ticular data set, <5c = 0.06, so that ~ w^, with ~ ~ L±/dv = 1.11 x 10^ s. On the 
other hand, ?adv = (1 AU)/?7 = 4.26 x 10^ s, which marginally satisfies the condition ?adv ^ '^kf 
However, if w^/vv^ were instead ~ 3 at 1 AU, as in the results of Bavassano et al (2000), then 
would be somewhat larger than ?adv Thus, the inhomogeneity of the solar wind may influ- 
ence the effects of cross helicity on the spectral indices. However, more detailed modeling of 
inhomogeneous solar- wind turbulence is beyond the scope of this paper. 

9. Comparison to Other Studies 

In this section, the results of this paper are compared to two recent studies of strong anisotropic 
incompressible MHD turbulence with cross helicity. 



9.1. Lithwick, Goldreich, & Sridhar (2007) 

The model for the cascade of energy to larger ^ used in this paper is based on the results of 
Lithwick, Goldreich, & Sridhar (2007) (hereafter LGS07). As a result, in both studies, if w+ and 
have comparable correlation lengths in the direction of Bq at the outer scale, then ~ ^ at 
all smaller scales. The principal difference between this paper and LGS07 lies in our treatment of 
the cascade time for the dominant fluctuation type, w+. LGS07 argue that if » w^, ~ 1, 
and <^ 1, then the shearing applied by w~ wave packets on a w+ wave packet at perpendicular 

scale kj^ is coherent over a time {k±w'^^)''^, which greatly exceeds the time (^|| va)~^ required 
for a w+ and w~ wave packet at perpendicular scale kj^ to pass through each other. In contrast, 
in this paper, it is argued that the coherence time for the straining of the w+ wave packet is of 
order the "crossing time" (^|| va)~^- As a result, the results obtained in this paper for the inertial- 
range power spectra, degree of anisotropy, cascade time, and energy fluxes are different from those 
ofLGS07. 

The approach taken in this paper is motivated by the following argument. As argued by 
LGS07 and Maron & Goldreich (2001), the w~ fluctuations propagate approximately along the 
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hypothetical magnetic field lines obtained from the sum of Bq and the magnetic field of the w"*" 
fluctuations. Let us call these hypothetical magnetic field lines the "w"*" field lines," and let us 
consider a wave packet of perpendicular scale kj^ and parallel scale (k\\)^\ where k\\ = ~ 

. Let us work in a frame of reference that moves at speed va in the —z direction along with the 
fluctuations. Let us also take an initial snapshot of the turbulence at t = and trace out all of the 
"w^ field lines" that pass through our wave packet. The volume filled by these w+ field lines is the 
"source region" from which the wavepackets encountered by our w+ wave packet originate. If 
we wait one crossing time {^va)~^ and take a new snapshot of the turbulence, then at any given 
location the w+ fluctuations will not have changed very much, since <^ . However, if we 
trace out the new w+ field lines passing through our wave packet, the volume that is filled by these 
new w+ field lines will differ substantially from the initial source region at distances > (^||)~^ 
from our wave packet due to the rapid divergence of neighboring field lines in MHD turbulence. 
In other words, small local changes in w"*" lead to large changes in the connectivity of the w"*" field 
lines. 

To see this, let the w"*" field line that passes through some point P in our wave packet at ? = 
be called "field line A." Let the field line that passes through point P att = (^||Va)^^ be called 
"field line B." As before, let us work in a frame of reference that moves at speed va in the —z 
direction. Field lines A and B are fixed curves, since they are traced out within two snapshots of 
the turbulence. If we follow field line B for a distance <^ it will separate from field line A 

by some small distance x that is 'tik^^. If we continue to follow field line B, its separation from 
field line A is analogous to the separation of two neighboring field lines within a single snapshot 
of the turbulence. As shown by Narayan & Medvedev (2001), Chandran & Maron (2004), and 
Maron, Chandran, & Blackman (2004), if a pair of field lines is separated by a distance x <^kj^ 
at one location, then the distance the field-line pair must be followed before it separates by a 
distance kj^ is a few times the parallel size of an eddy of perpendicular size kj^ - i.e., a few 
times (^11 )^^. It turns out that the particular value of x has little effect unless one considers the 
(irrelevant) case in which jc/J ~ e^^, where A' is large and d is the perpendicular dissipation scale. 
(Chandran & Maron 2004) This is because within the inertial range the amount of magnetic shear 
increases towards small scales; therefore, if x is made very small, then the distance one has to 
follow the field-line pair in order for x to double becomes very small. Thus, as a result of the rapid 
divergence of neighboring field lines in MHD turbulence, the "source region" of our w"*" wave 
packet att = {^va)~^ differs substantially from the source region at ? = at distances > (^||)~^ 
from our wave packet. Because the w~ fluctuations vary rapidly in the direction perpendicular to 
the magnetic field, the w~ wave packets encountered by our wave packet will decorrelate on a 
time scale of order the crossing time {^va)~^ due to the time evolution of the source region. 

It should noted that there are two unexplained aspects of LGSOV's model, as pointed out 
by Beresnyak & Lazarian (2007). The first concerns the nature of the transition from the weak- 
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turbulence regime <^ 1 and ^ 1) to the strong-turbulence regime ~ 1 and ^ 1). 
In LGSOV's analysis, as one passes from the weak regime to the strong regime, the coherence time 
for the straining of the w+ wave packets by wave packets increases by a factor of (%^)^^ and 
the energy cascade time decreases by a factor of . It is not clear why these large changes 

—5 /3 

should occur across the transition scale. The second issue is that LGS07 find that E'^{k^) oc 

_c /3 

and E^{k^) oc k^ regardless of the fractional cross helicity. Since the ratio E^{k^_)/E^{kj_) 
is independent of wavenumber, it is not clear how pinning could occur in their model, or how the 
spectra would behave near the dissipation scale. 



9.2. Beresnyak & Lazarian (2007) 

Beresnyak & Lazarian (2007) (hereafter BL07) have published an online article on strong 
MHD turbulence with cross helicity. The following discussion refers to the version of their article 
that is available electronically as of the writing of this paper. The work of BL07 is similar to this 
paper in that both studies take the dominant fluctuation type, w~^, to undergo a weak cascade. Also, 
equations (|3T1) through (|33l) of this paper are equivalent to their equation (5), except for the fact that 
they take the parallel correlation lengths of w+ and to differ by a constant multiplicative factor 
when a power-law solution for the spectra is assumed. (The possibility of more general solutions 
is claimed by BL07.) On the other hand, there are a number of significant differences between this 
paper and Beresnyak & Lazarian's (2007) work. They argue that for the w+ fluctuations, the dom- 
inant nonlinear interactions are between fluctuations with comparable parallel correlation lengths 
and different perpendicular correlation lengths, whereas for the w fluctuations the dominant in- 
teractions are between fluctuations with comparable perpendicular scales. Here, it is argued that 
for both w+ and the dominant interactions are between fluctuations with similar perpendicular 
scales. When ^ w^^ and ~ 1, their Figure 1 suggests that the parallel correlation length 
of fluctuations can be less than the parallel correlation length of w~ fluctuations. It is argued in 
section |7] of this paper that this can not be the case, because the w+ fluctuations will imprint their 
parallel correlation length onto the fluctuations. They argue that the scalings given by equa- 
tions (|3T]) and (|33l) [equivalently, their equation (5)] can not apply if the parallel correlation lengths 
of the w+ and fluctuations are equal, arguing that this would require £+ = whereas in this 
paper the ratio £+/£^ depends upon the slopes of the power spectra, as in weak turbulence. They 
argue that if the w+ and fluctuations are driven with the same parallel correlation length at the 
outer scale, there will be a non-power law part of the solution at large scales that will transition at 
smaller scales to a power law solution with and both oc kj_^^^ and k^ oc k^^ . In contrast, 
in this paper a power-law solution starting at the outer scale is obtained with different scalings for 

and , and with growing more slowly than k^ . 
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10. Conclusion 

This paper proposes a new phenomenology for strong, anisotropic, incompressible MHD tur- 
bulence with cross helicity and introduces a nonlinear advection-diffusion equation [equation (fT5l) 1 
to describe the time evolution of the anisotropic power spectra of the w+ and w fluctuations. It 
is found that in steady state the one-dimensional power spectra of the energetically dominant 
fluctuations, E^{k±^), is steeper than , and that E'^{k±^) becomes increasingly steep as the 
fractional cross helicity Oc increases. Increasing Oc also increases the energy cascade time of the 
w+ fluctuations, reduces the turbulent heating power for a fixed fluctuation energy, and increases 
the anisotropy of the fluctuations at small scales. 

Although most of the discussion has focused on forced, steady-state turbulence, the results of 
this paper can also be applied to decaying turbulence. For example, equations (fT2l) and (fT3l) can 
be used to estimate the time scale for turbulence to decay. The resulting prediction is that if the 
fluctuations are initially excited with ^ and with comparable parallel correlation lengths 
at the outer scale, then the turbulence will decay into a state in which fluctuations are absent, as 
in the earlier work of Dobrowolny, Mangeney, & Veltri (1980), Grappin et al (1983), and Lithwick 
& Goldreich (2003). This "maximally aligned" state will then be free from nonlinear interactions, 
and will persist for long times until it damps via linear dissipation. 
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